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ABSTRACT 


This thesis investigates the performance degradation resulting from 
multitone interference on orthogonal, noncoherent, frequency-hopped, M-ary 
frequency-shift keyed receivers (FH/MFSK) where the effect of thermal and other 
wideband noise is not neglected. The multiple, equal power jamming tones are 
assumed to correspond to some or all of the possible FH M-ary orthogonal 
signaling tones. Furthermore, the channel is modeled as a Ricean fading channel, 
a possibility precluded when thermal noise is neglected. Both the signaling tones 
and the multiple interference tones are assumed to be affected by channel fading. 
Both band and independent band multitone interference are considered. 
Performance is evaluated by obtaining a union bound on the probability of bit 
error, and receiver performance is compared with exact results for band multitone 
interference of a noncoherent FH/MFSK receiver under comparable 
circumstances. Except for the case of Rayleigh fading of the signal, the union 
bound is very tight for those cases that can be compared with exact results. The 
advantages of the union bound approach are twofold. First, the union bound 
approach yields a solution that is far less computationally intensive than that 
obtained with the exact approach. Second, the union bound approach allows 
numerical results to be obtained for interference conditions that are not amenable 


to exact analysis, such as independent multitone interference of FH/MFSK. 
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I. INTRODUCTION 


A. SCOPE OF THESIS 


The goal of this thesis is to determine the electronic counter-counter measures 
(ECCM) potential of a non-coherent frequency-hopped M-ary frequency-shift keying 
(FH/MFSK) communication system under conditions of worst case hostile electronic 
counter measures (ECM) and fading channels. The performance of the FH/MFSK 
receiver will be investigated for frequency-nonselective, slowly fading Ricean channels 
and hostile multitone interference. 

The receiver considered uses a quadratic detector in each M-ary branch of the 
receiver for signal detection. Both band and independent multitone interference are 
considered and both jammer and signal are assumed to be affected by the fading channel, 
where it is also assumed that the channel fading need not necessarily affect the signaling 
tones and the interference tones in the same way. The cases of Rayleigh, Ricean, and no 
fading are all considered. 

Performance is evaluated by obtaining a union bound on the probability of bit 
error. The results are compared with the exact results, examined in [1], to determine the 
accuracy of the union bound approach. Then the union bound is used to obtain probability 
of bit error for situations where the exact approach is either very difficult or even 


impossible to obtain. 


B. IMPORTANCE OF A SPREAD SPECTRUM COMMUNICATION SYSTEM 


The benefit from investigating the performance of a spread spectrum 
communication system lies in the extent to which spread spectrum communications are 


currently used. Their application to any military communication system aids in achieving 


secure, low probability of intercept, and/or low probability of detection information 
exchange between friendly units as well as reducing the ability of hostile elements to jam 
friendly communications. Other applications of spread spectrum are: 

1) The use of either frequency-hopping or direct sequence spread spectrum in 
multiple access communications (as in cellular phones where vehicles communicate with 
a central station) in which a number of independent users are required to share a common 
channel without an external synchronization mechanism. 

2) The application of spread spectrum techniques in wireless local area networks 
(LANs). Wireless LANs are very important for the establishment of a network for mobile 
forces that circulates among the cooperating units various kinds of information, such as 
text messages, video, or audio. Wireless LANs are also very popular in commercial 
communications and are used for wireless communication of laptop computers with a 
control station. 

In any case, the potential of frequency-hopping spread spectrum to reject 
interference of either unintended or intended transmissions with the communicating 
signal is of crucial importance. The problem considered in this thesis is the evaluation of 
the performance of a frequency-hopped system in order to determine the system's 
capabilities and limitations. 

In chapters 2, 3, and 4, the union bound approach to performance analysis, a brief 
description of frequency-hopping, and the notation used in the thesis are discussed. In 
chapters 5 and 6, receiver performance with independent and band multitone jamming, 
respectively, are analyzed. In chapter 7 numerical results for both independent and band 


multitone jamming are presented. 


Il. IMPORTANCE OF UNION BOUND APPROACH 


The union bound approach simplifies the mathematical development by deriving 
an expression for the upper bound on the probability of symbol error that is always much 
simpler than the exact form. There are usually two approaches used to evaluate the 
probability of symbol error of a communications system. 

In the first approach, known as the indirect approach, the probability of symbol 


error is obtained using the relation 
aes (1) 


where P. is the probability of not making a symbol error. Noting by X,,X,,...X,, the 


random variables representing the output of each branch of the receiver for a signal with 


M symbols, we have: 
Pee ey GK, > Xf)... K, > (2) 


The calculation of this probability is difficult, especially in the case where we deal not 
only with gaussian noise but also with jammer tones affecting the signal. This will 
become particularly obvious later on when we analyze receiver performance with 
independent multitone interference. In that case, the joint probability distribution 
functions required to solve the problem using the exact solution approach make it 
impossible to obtained a closed form expression for the probability of symbol error for a 
system using more than two symbols. 

However, if we use the union bound, or direct approach, the solution works as 


follows, providing an upper bound to the probability of error calculation: 


see = !27 2G Oy , -- OEe < Xyll) S Pr(X, <X,[F Pe, i, 
+...+Pr(X, =< X 41) 


We see that in place of joint probability distribution functions (pdfs) we now have the 
sum of individual pdfs. In fact, the individual pdfs can be grouped in groups of identical 
pdfs occurring with a certain multiplicity, so that the solution is a linear combination of a 
few pdfs. 

Summarizing, the union bound approach provides a simple, fast method to 
investigate communication and jammer signal scenarios that are very difficult, if not 
impossible, to explore with the exact result approach as described in [1]. By using the 
union bound, we can increase our knowledge of communication systems performance for 
a larger variety of jamming conditions. However, since this evaluation is a bound, its 
accuracy, or tightness, needs to be tested so that we can comment on its effectiveness and 
limitations as a tool. In fact, in the following pages it will be shown that in some cases the 
union bound approach is so accurate that it provides virtually identical results with the 
exact solution, while in some other cases the results obtained display an anomalous 


behavior. 


Iii. BACKGROUND 


The main idea of frequency-hopping spread spectrum communications is to 
provide protection against a hostile jammer by increasing the bandwidth that the 
information signal occupies much more than the minimum required. By increasing the 
frequency range over which a system effectively operates, we force the jammer to spread 
its power over a wider frequency band and thus be less effective from a "per-frequency” 
signal corruption capability point of view. 

The main idea of frequency-hopping spread spectrum is based on the 
multiplication of the signal prior to transmission by an intermediate frequency, generated 


by a frequency synthesizer, of the form 


f, =f, +(i-1)-Af,, where i=1,2,...N (4) 


where N is the maximum number of possible frequency hop bins, Af, is the separation 
between the carrier frequencies of adjacent bins, and '1' changes pseudo-randomly every 
T, seconds. By so doing, the entire spectrum of the signal transmitted, is shifted from its 


carrier frequency f, to the new carrier frequencies: 


eet. . (5) 


The result is that a new data modulated carrier that is shifted from one frequency to the 
next 1s created. From (5) we see that there are N different frequency hop bins, each of 
bandwidth Af, . The value of i (and thus fc; ) is changed periodically according to some 
predetermined, but apparently random to a third party observer, noiselike spreading code 


called either the pseudo-random or the pseudo-noise (PN) sequence. 


Generally, frequency-hopping is used with noncoherent modulation formats such 


as noncoherent MFSK. Such a signal is described by: 
s(t) = /2A, cos|2n(f, +(m—1)-A,)t+ @ | (6) 


where m=1,2,...M, 9, is the i" symbol phase, and A, is the frequency separation between 
each of the M signaling tones. Multiplying this signal by a signal at the intermediate 


frequency f, = f, +(i-1)- Af, , we get: 


s'(t) =2/2A, cos{2n]f, +(m-1)A,|t+9, \ cos(2nf,t) = 
SAPS, cos} 2a f; +f. +(m- 1A, |t +0, ; +V2A. cos} 2a f, -—f.-(m—-1)A, It +0. \ 


(7) 


The carrier frequencies of the first part are 
f,+f,+(m-l)A, =f, +(i-DAf, +f +(m-I)A, (8) 


and is smallest for i=] and m=]; 1.e., f, + f,. The carrier frequencies of the second part 


are 
f, -f, -(m-1)A, =f, +(i- Af, -—f, -(m—-1)Af (9) 
and is largest for i=N and m=]; 1.e., 


f +(NEDAGe f.. (10) 


If the smallest frequency of the first part is greater than the largest frequency of the 
second part, then 


2f, > 2f,+2B+(N-l)Af, => f, > f, +B+(N—-I)Af,, /2 (11) 


where B is the required guardband above and below the high and the low frequency 
signaling tones respectively. In this case, the signal s’(t) can be high-pass filtered to 


remove the frequency difference contribution, and the frequency-hopped signal becomes: 
V2A, cos{2n[f, +f, +(m-1)A,|t+9,} (12) 


One obvious advantage is that a hostile jammer needs to know the pseudo-noise 
sequence of the communicator's transmitter in order to jam the specific hop bin where the 
transmitter operates at each time instant. To overcome a follower jammer, the signal must 
hop to a new carrier frequency before the jammer is able to complete its tracking process. 

Another advantage is due to the fact that the power spectral density of the 
frequency-hopped signal is identical to the that of the conventional signal in a specific 
hop bin. However, since the signal hops from bin to bin, and assuming that the 
probability that any bin is occupied is equal to (1/N), the average power spectral density 


1S: 


Mz 


PSD= 


Sl Se Ie (13) 


ie 


. 


Ik 
—- 


Hence, the signal power spectral density is lowered by a factor of N on average. 
Before proceeding we need to distinguish between the two main tone jammer 
methods used against communication systems. The first is band multitone interference 


and the second is independent multitone interference. In each case, we assume that the 


multiple jamming tones are all of equal power. In band multitone interference, the 
multiple equal power jamming tones are distributed randomly across the entire frequency 
hop bandwidth, but the number of jamming tones placed within each jammed frequency 
hopped band is specific. In independent multitone interference, the multiple equal power 
jamming tones are distributed in a completely random fashion, meaning that the number 
of jammer tones within a frequency hop bin can vary from zero to a maximum equal to 
the total number of signal tones corresponding to each symbol the communicator can use. 
We also assume that each jamming tone corresponds exactly to one of the possible 
M x N signaling tones. From the communicator's viewpoint, this is the worst case. 
Before proceeding to the analysis of the receiver performance with multitone 


jamming, the notation used in the remainder of the thesis 1s discussed in the next chapter. 


IV. FH/MFSK SYSTEM, CHANNEL MODEL, AND SYMBOLOGY USED 


The FH/MFSK communication system examined is assumed to have N non- 
overlapping FH bins. The bandwidth of each one is B Hz, where B is the bandwidth of 
the non-hopped MFSK signal. So the total frequency-hopped bandwidth is augmented to 


be greater than or equal to B-N. We also express by M the modulation order where 
Miao (14) 


and k is the number of bits per M-ary symbol. Since the number of signal tones within a 
hop bin is M and the total number of bins is N, the result is a total of M-N possible signal 
tone locations. 


The symbol duration T, is related to the hop duration T, by: 


: (15) 


where K=1,2,3... . The symbol and hop rates are the reciprocal of the corresponding 


durations, so 


pot 
} 


S 


(16) 


l 
and R. =—., i 
ac; (17) 


respectively. Symbol rate is related to hop rate by 


R,=K-R, (18) 
Using the frequency hopped signal (12), we note that the M symbols are assumed 
to be equally likely, both f; and f, are taken to be integer multiples of the symbol rate Rg, 
and the frequency difference between signal tones Aris taken to be an integer multiple of 


the symbol rate. Hence, 
A, =p:R, (p integer). (19) 


To incorporate the idea of fading channels, we model the amplitude A, as a Ricean 
random variable. This includes Rayleigh fading as a specific case. 

The transmitter and receiver are assumed to function in perfect synchronization, 
and the receiver after dehopping the signal uses a bank of M quadrature detectors where 
the integrator time constants are normalized to the symbol duration for notational 
convenience (Figure 1). 


The existing noise interference is expressed as AWGN with a flat average power 


spectral density defined as x . The total power of the jammer is expressed as P,_, and 


the jammer delivers q equal power interfering tones such that the power of a single 


interfering tone is equal to 


P 
P.=—. (20) 





The subscripts used in the description of the Ricean random variables are ‘c' for 


the signal parameters and 'j' for the jammer. Thus, the pdf of A, is 





2 Z 2 2 
f, (a,) = - on ete) 1f 222] -u(a,) (21) 


2 
VOwe 


10 


where u(¢) is the unit step function, I,(¢) is the modified Bessel function of the first kind 


and order zero, a*. is the power of the direct signal component, and 207. is the power of 
the diffuse signal component of the respective tone. The total average power of either a 
Ricean faded jammer or signal tone is equal to the mean squared value of the amplitude a, 


or a;, respectively, and is 
P. =a’; +207; (i=c,j). (22) 


The average power is assumed to remain constant from hop to hop. It is convenient to 


define the signal-to-noise power ratios of the direct signal components as 


2. 


p, = 24 (i=c,j) (23) 
Oo 


n 





and the diffuse signal components as 





B= (isa) (24) 


where 6”, is the noise power. 
For orthogonal signaling, the average probability of bit error is related to the 


average probability of symbol error by [2] 


M /2 
es oesliea 25 
i (25) 


1] 


and the average energy per bit is related to the average energy per symbol by 


E, = log,(M)-E, =k-E,. (26) 


V. FH/MFSK RECEIVER PERFORMANCE WITH INDEPENDENT 
MULTITONE INTERFERENCE 


With independent multitone interference, the multiple, equal power jamming 
tones are distributed randomly across the entire frequency-hopped bandwidth. Two 
possible cases exist concerning the signal and jammer tones within a specific frequency- 
hop bin. In the first case there is no interference at all in the bin (the jammer tones are 
elsewhere), and the probability of error is the same as for a conventional MFSK signal in 
the presence of AWGN only. In the second case there is jammer interference in the bin. 
Since a bin is composed of M different signal locations and a potential jammer may 
place one or more of its tones to the aforementioned signal locations, this case is further 
divided into two subcases. In the first subcase the signal and jammer tones are on 
different branches, while in the second subcase one jammer tone coincides with a signal 


tone. This gives a total of three individual cases that are treated below. 


Case |: 

In the following, we express by R,, R,,...R,, the random variables representing the outputs 
of branches 1,2,...M, respectively, of an MFSK receiver. Without loss of generality, we 
assume that the symbol '1' is transmitted. In order not to have an error, the output of 
branch one must be greater than the outputs of all other branches. The opposite situation 
expresses the error condition. As in equation (2), the probability of symbol error provided 


by the union bound is 
P, = Pr(R, <R, UR, <R,U...UR, <Ryyll) 


<Pr(R, <R,|l)+Pr(R, <R,ll)+...+ Pr(R, <Ryll) 


(27) 
<(M-1)-Pr(R, <R,|l) 


and the last equation is true since the individual probabilities are equal. 


Case 2a: 

Suppose that we have L jammer tones in the hop bin where the signal is. By S, we express 
the existence of a jammer tone, and the subscript after ] corresponds to the specific branch 
where the jammer tone is present. For example 8; means that a jammer tone is present in 
branch one. To simplify the discussion for this case, suppose that we have two jammer 


tones (L=2). Now 


Pr(error|1,S, ,S;)< Pr(R, < R,/1,8,)+ Pr(R, <R,|1,S; )+ Pr(R, < R,|I)+...+Pr(R, < Ryll) 
Since Pr(R, < R,|l) = Pr(R, <R,|l) =...=Pr(R, < Ry |D= 


= Pr(error with no interference in signal branch) (28) 
and Pr(R, <R,|L,S; ) = Pr(R, < R,|LS, ). (29) 
then (28) reduces to 
Pr(error|l,S; ,S;)<2Pr(R, < R,|LS8,,)+(M—3)-Pr(R, < R,|1) (30) 


Generalizing for L jamming tones per hop bin, 1<L<M-1, we get 


IE 
Pr(error|l, )1S,_) < L-Pr(R, <R,|L,S, )+(M-L-1)-Pr(R, <Ryll) (31) 
j=} 


Case 2b: 
As with case 2a, we begin by considering two jammer tones. In this case 
Pr(error|l,S; ,S,,) < Pr(R, < R,|LS8; .8,)+Pr(R, < R,I1,S; )+Pr(R, < R,ILS, )+-.. 
+ Pr(R, < Ry ILS; ) 
(32) 


As before, this can be simplified to 
Pretror|l;s, .S, )= Pr(R, <R,/1,S, ,S, )+(M@—2)- PréR, < RAILS, ). (33) 
Generalizing, for L jamming tones per hop bin, 1<L<M-1, we get 


L 
Pr(error|l,S, »s) =(L-1)-Pr(R, < R15, .5,)* (M-L).- Pr(R, = RwllS, } 
]=2 


where it is understood that the summation is zero for L = 1. 


(34) 


Now the probability of symbol error with at least one jamming tone per hop bin (L21) is 


given by: 


L L L L 
Pr(error)= Pr(error]l, » eels » S,_,)+Pr(error|l,S, , > SF os yds 55.) 
l=] j= l=2 2 


(35) 


Since we have L jammer tones and M signaling tones per jammed frequency-hop bin, 


L 
Pr(1,S,,/S, == (36) 
] = 1+] M 
iC 
and Pr(1,>_S, eonepmeliaa (37) 
lel I+! M 


So the probability of symbol error for L jamming tones per hop bin, L21, is given by 


IL 
P(L)< — (L -1)-Pr(R, <RjILS, ,S,, ]+(M-L)-Pr{R, <RylLS, )|+ 
IL. 
- aa; JL-Pr{R, < R,|LS,, }+(M -L- 1)-Pr(R, < Rll) 
(38) 

We now look for the probability of having L jamming tones per hop bin. This 
leads to the expression for the total probability of bit error. For independent multitone 
interference, and since we assume q jamming tones randomly interfere with a total of 


N-M possible signal tones, 1<q<N-M, the probability that a frequency-hop bin contains 


one interfering tone is 


ie a 
nave (39) 


The probability that a frequency hop bin contains a second interfering tone 1s 


sj 
— * (40) 


The probability that a frequency hop does not contain an interfering tone is 
q 
a 41 
<i (41) 


The probability that a frequency hop bin does not contain a second interference tone is 


i (42) 


16 


From (39) - (42), we derive the probability that a frequency hop bin contains no 


interference tones as 


ror (stab) (-wita)-C-eataca 


(44) 


In the same manner, we define the probability that there are two jamming tones in a 


frequency hop bin as 


yay waa) wha) 


-Ugaa -. (- M- axa. 


i=] 


(45) 


Combining (43), (44), and (45), we get the probability that L jamming tones are in a hop 


bin, where 1<L<M-1, as 


17 


Pr(L jamming tones per hop bin) =| .. 
q-(L-1)-1 q-(L-1)-1 
M-N-(L-1)-1]"" | M-N-(M-1) 


[esa] Scs) 
=, MeN es M-N-i 


(46) 


The situation where L=M must be treated separately. Looking at the two 
individual products multiplied together in (46), we see that the first product has L terms, 
corresponding to the probability that L out of the M signals of the bin are jammed, while 
the second product has the remaining M-L signals of the bin which are not jammed. Since 
when L=M all the signals are jammed, for that case we need only the first product term. 


So the probability of having L=M jamming tones in a hop bin is 


L-1 


q-1 
Pe —— 
see So vaer 


(47) 
Now, the total probability of symbol error for 0<L<M is given by 
Pr(error) = Pr(L = 0)-P,(1) + Pr(L = 1)-P,(1|L = 1)+...+ Pr(L = M)-P,(1|L = M) 
(48) 
The maximum number of jamming tones we can have in a bin is M if M<q 


and q if q <M. This can be represented as min(q,M). By combining equations (46), 
(47), and (48), the probability of symbol error is given by the relation 


1= 


{Ets st] 


P.< Ht: -s)} ~1)Pr(R, <R,|l) + 


L~1 


q— 


step 209 


(49) 


The probability P,(M) is obtained by letting L=M in (38). This yields 

P,(M) =(M-1)-Pr(R, <R,|1,8; ,S, ). Based on the above probabilistic description, it 
is apparent that we need to obtain the following probabilities for the MFSK receiver: 
Pr(R, <R,/1), Pr(R, < R2/1,S;, ), Pr(R, < RAILS, ), 

and Pr(R, < R,|1,8, ,8;, )- 


A. Pr(R, <R,|1) 


For noncoherent, orthogonal BFSK, Pr(R, <e 1) is equivalent to the probability of 








bit error: 
1 
Pa.) = 7 aa (50) 
N 
where o”, = ae (51) 
b 


For noncoherent, orthogonal MFSK, the only thing that changes 1s that T, is replaced by 
T ° 


Ss 


N 
a 2 
o 1 (52) 


To account for the existence of a Ricean fading channel! we multiply equation (50) with 


(21), the Ricean pdf, and then integrate over all possible values of a,. This yields 





“1 Gaius ae tO, a hat 
P= iP (a.)-f, (a,)da, = I> Pe | = 26?, i em aia ye a 5 -u(a, ) 


Cc 


- ] ~a,2 i | lan l 1o( 2} 53 
oe eee oe ee e°) 


Let x=a*, => dx = 2a,da,, and (53) becomes: 


] a” | 4 ] avx 
r= geen if fs em|-§[ Soo yf Ep la a.. (54) 
Cc Cc 0 Nn c Cc 


We then use the identity [3]: 











ee m!p” = | : 
Jx 2 .e J, (2BVx )dx = Le 


(55) 


where 3%",,(e) is the Laguerre polynomial and is defined [3]: 
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BY 
ey" a 
In our case m=n=0 and the whole polynomial evaluates to unity. In place of "a" we have 
1 + Lo 
Nok o* 
] oe) 
1 Em 43 ms) Jexp(- ~ax)-J »(2BVx )dx = 
| ie | ( Peet | ex les ;)| 3s83 
a ie we eXp i ee aE at eee — 
40°, Oo. +0 40“ )\o,° +0°, 


Oo’, = l cme (57 
li 2, exp IG eel ) 
2+2:- ae 


wo in place of "B" we hav 





? 





Cc 














ee 2 
O.+O 


n 


Using (23) and (24), we have the definitions: 


é ; a 
p.= direct signal-to-noise-ratio = — 





n 











20° 
€. = diffuse signal-to-noise-ratio = = 
Applying these definitions to (57), we get 
| P. : 
Pr(R, <R,|l)= ee ex] 52. (58) 


a 


B.Pr(R, <R,{1,Sj,) 


This probability of error is the same as that derived in section 4.A except that a 
jammer tone is on the signal branch. Since signal and jammer tone exist on same branch, 
we can consider that they constitute a new kind of (altered) signal existing on the branch. 


So the probabilistic analysis yields the same result with the only difference being that the 


"old" signal power a’. is replaced by 
a*, +a°; +2a,a,cos0,. (59) 


In other words, the power of the signal and jammer tones add together but not 


necessarily in phase. We then integrate with respect to 0; after multiplying by its pdf to 
eliminate the dependence on 9. The phase angle®, is assumed to be uniformly distributed 
in [0,27] ; 

Without jamming, the mean values of the random variables representing the outputs of 


the in-phase and quadrature components of the first branch are given by: 


X= V2a, cos(6) (60) 
X,, = ¥2a, sin(6) (61) 


and in this case the output of the first branch is given by: 


X?,, + X74, = 2a”.(cos?(6) + sin?(@)) = 2a”. (62) 


This results in the probability of error: 
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l aac. 
Pr(R, <R,|l) = 5 exp crema (63) 


With jamming we have for the in-phase, quadrature components, and their squared sum 


as the final branch output the following relations, respectively, 


X,, = V2a, cos(9) + v2a, cos(9' (64) 
— = V2a, sin(@) + V2a, sin( 0} (65) 
and X?, + X74, = 2a”. +a”, +2a,a, co 0") (66) 


We define 0, = 0 — 0; as the phase of the jammer with respect to the received signal to 


get the final form of equation (59). Considering the equations above, we can understand 


why we substitute the previous signal power a’, in Pr(R, <i 1, Sj, by the composite 


power a’, +a; +2a.a cos{@ i). Substituting (59) into (63), multiplying by 
l ’ . 

fy, (6 ; =e and integrating over 0 ,, we get 

a 

4n 


2n - 2. 2 a bE 
Pr(R, <R,|LSj, ) = joaf eee, (67) 


2N, 


We then use the identity 


2n 


i) =~ Jexp|+ Z: cos(o - 6) kde (68) 


to obtain 








c 


‘ | Ae a’; ace aj 
Pr(R, <Rall,Si,) = 5 exp] - 55 1 epg 





] ] aad 
= Fe] sa le +a’) | 9 | (69) 


On 


However, this probability of error is conditional on both a, and a,. To eliminate the 


dependence on a,, we multiply by its pdf (Ricean) and then integrate for all possible 
values of a,. After multiplying (69) by (21) and integrating, we get: 


Pr(R, < R, ILS), = 


l l i l l ae +0". a.°a; a. -O, 
pee aac tl) a en[-aetedf eof Sse ]C32) Se 














l ae on 
= FoF, eHP|— Z57, (27s) |-expl— 9 | 








(70) 


We now have the identity[3]: 
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. ] Le 
fx-evpl-ax"]-1,(0x)-4,Crn)ox= 3 eng = ay (71) 


The correspondence of the terms in (70) with the ones in the identity are: 


72 A 
Oo +6'n g a, jo 
= 2 =o > Se" > - 
20°.07, ” Co, on 





a 





So the probability of error, after eliminating the dependence on a,, is 


Pr(R, <R BS.) =e ty eng -3 S| to 0), (72) 
Sede 2(c?, +07, ) 2 0. +07, oto, 2 


To eliminate the dependence on a;, we multiply by its pdf (Ricean) and then integrate for 
all possible values of a;. The Ricean pdf for a; is given by: 


; | ) a, a ta,” Oj . ( ; 7 








2 


2 2 2 2,2 P] Z 2 2 
5 . 


(67. +o°,) oy = (oer +67, )o?, 


( om a,}-1 aes 
j 0 2 
“\a’.t+a*, oe 





] 
2 02,(o%, +07) oP = (c? Gs Aloe, 


Z Pi 2 2 
| 07,07. +(o jO°6 +A jO . ( O, oeuae 
0 SS | eee 


oa a; l oe ) 





(74) 


To integrate (74) with respect to a, we again use identity (71). The correspondence of the 


terms in (74) with the equivalent ones of the identity are: 


57, +(67, +07,) a mc 
A el \>)|(S = yy ee oe 
z (oa los OctOn Oj 


For the argument of the Bessel function, we have: 


OL | “4 
2 en (ee aa 
By Oe t+On VO’; OO, 


iz = Jie 2 2 
2a Ab 57, +(o7, +07.) O°; +O°.+0'n 
2 2. \eeee 
2 (o?. +6 »)o?; 


C2 
So the final Bessel function has the form: I {24 (75) 
j+o .+90 
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So the terms appearing in front of the integral of (74) reduce to: 


l 1 ie in 
(.+04)-04,| aig tyeeties (76) 


a 2 2 rs = 
of 2 3?(o7, +07.) IN \4 hoy +07n} 


Finally, the argument of the exp(e) function that the integral of (74) yields is: 


(85) (i 4) 
B?-y? or onan J GS 
4a - E “eset 


2 AG? toils’; 





and combining this with the leading exponential functions we get the final exponential 


function to be of the form: 


2 
= (i .| 
j 20S? 
Oe t+O'n Oey ta jo, +(a joe oul 


] 
| 5 jeg ee = (g*, +67, )o7; = iG eo ae (77) 
De lot be i 


Combining (75), (76), and (77) as required by (71) we get: 


si) 


ya Z 
oer a7. +07; 


Pr(R, <R,[1,Sj,) = 











Substituting 
oO. a : Go : 20°; 
P. o.° P; mm Gt? 7g, 


in (78), we get 


] P. TP; Z PePj 
Pr(R, < R,|L,Sj,] = ey, =e E+E, i) ule +. 


C. Pr(R, <R,|LS, ) 


eT a a }al 
-ex —_.S |. 
2(6?, +07. +07) a 2 O°; +0°.+On,/ ° 


OO, 


OSG Sb 


(78) 


(79) 


We begin by designating V, and V, to be the random variables representing the 


outputs of branch 1 (where the signal resides) and branch 2 (where the jammer tone 


resides), respectively. Their pdfs are given by 








] v, +2a%¢ Jol 2v, 
f [v ll,a |- 2X) — See | eee 
WV eed aa 2 
267, 2674 o- 1 
1 V+ 2a j a. 2V, 
i, (v,lla,) = 7 XP] — 5 ‘I, 
2 20 n 20"n On 
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(80) 


(81) 


| 


Multiplying (80) by (21) and integrating, we get 


albdlemtaredets) ae 
V, V1 ~ 907,67, PL 267.) *PL 26? 


if exo a’ Ta the.) I aa 2V, 1 2a 
Jee p Sy 0| G2 "Io ao a. 























(82) 
We then use the identity given in (71) where 
Gan 1D) Sere ZY. a, 
— 5. B= ot La (85) 
Now we have 
nia Coe 
Onset oe 2v,\Q, 
TE os", Veo Ol (84) 
C1” 
Zr Ol 
B-y? of, a v0 + 0%, re 
4a eee o-.) 207, +20 .\o'.0°- — 
or ome 


Hence, using (83), (84), (85) in (71), we can evaluate (82) to obtain 


¢ ( 1) It ( l v, oe | ZY PeOhe 96 
V = OO _—_—. oo Ce = —— + 
oS (67, H 27.) we 2 o', +200) “Won 22ore od 
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Similarly, 


6, (vals, ! [4 e201) | Ata ah 
MAN} (G7, 4207;) Th 2 0%, +207)) %o*, +207, 
Using the random variable transformation R* = V,, we get 


I, lr F205. V2r, -., 
ee _ oo a eee 
fe, (Rill) (G7, +207.) ox 260,426.) Noma Ce 


and similarly, 


Jt I, ‘20 I ( V21, -o, | 
=e0 


I, 
fp, [R,IS,, = (c?, +207;) ox - 2 ; co’, +207; eee (89) 


O°, +207; 


Since R, and R, are independent random variables, the joint pdf of R, and R, is equal to 


the product of the individual pdfs: 


far, (r, oT, LS,, = fk, (r, 1): fR, (r, s,,] (90) 


a co 


So PR, <R,ILS,, }= | [fe, -fe,drndr, = [fp (8) (Si (91) 
0 


a] 


For the inner integral, 
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20 Fa: A 
ffs, (rolS;, Jar “lag on- site| 8 (92) 


2 2 
OG’, +207; Con 


2 2 
Ip dr, I, X-4fo7, +207; 
a = dk = and SO et—‘“‘i‘a dW 


et 2 2 > and —; a. i 2 
O°, +20°; 6°,+20 ; O°n +207; O°n +207; 


I, 


FS _—s—(94 
On +20°; os 


Also to transform the lower limit of integration, when r, =1,, x= 


Using (92), (93), and (94) to (91), we get 


20°; 
“ x LOweeeoau V2 +0, 
Pr(R, < Ry ILS, }= J X€XP| — >" — 5 sane x |dx 








(95) 
Now compare the integral above with the Marcum’s Q-function [2], [4]: 
¥ x? +a? 
Q(a,B) = ]x-exp| - ‘1, (ax)dx, (96) 
B 


and notice that if we let 


Oj Ir, 


3= (Sa ee 


3] 


Then (95) is given by 


aX =o EL a | 
ae occ, o7, +207; 


(97) 


Substituting (97) in (91), we get 


a4 2 2 
r laze. A 2a 
PR, <Rjll,S;,) = [pice om vn.) 
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98 
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Consider the following identity[5]: 





“1 a R R any mR aR 
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Goes al a Vo? ar) pete a7, 4+6072_ fate Vo? =| 
(99) 
If we let 


= 2 2. i ~S a 2 2 = ; 
6, =0',+20°;, a,=a;-v2, R,=4,, 6, = Vo',+20°, , a, =v2 Or 


then (99) can be used to evaluate (98) with the result 


op: 





a; 






Pr(R, <R, ILS, | Gato. |, Q a" 

r < . le .- = ee a 
! Aleman Ac?, +67, +075) o?, $07, +07; ; o?, aa teens 2 
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(100) 
Using (23) and (24), we get 
1+6. y2P. aes 
PRs < RIS, 2s oe I d Tb, 8p Pe | 
1+6; (20; V2 
eee Jae 8, 2+ +E, 
(101) 
We can now use the identity [5]: 
O(a.) =1- (6.2) +ex9{ -! 7 )-t6(0.) (102) 


in (101) to get 


Pr(R, <R,|1,$j,) = Jed Series 


1+é, sees P. TP; | I 2 PP; 
248,48, PU 248, +8,) 248, +6, 


(103) 


D. Pr(R, <R,ILS,_,S,, | 


This probability of error is obtained from Pr(R, <R,|1,S :) by altering the signal 
in branch one to include the jammer tone S 3, AS before, a*. is replaced with 


a’, +a’; +2a,a, cos( 0 a We must then integrate twice to eliminate the dependence on 


the fading parameters. The phase angle 8 requires a third integration. However, the 
integrations cannot be evaluated analytically. The triple numerical integration is time 
consuming. As will be explained later, an effort is made to approximate 9 so as to avoid 


the third integration. 


The rR, <KIES re | for a non-fading channel is given in [6] as 


Pr(R, < R,[1, Sj, ) a ah- - ay Eanes (104) 


Using (102), we get 


ja’ ja’, ja la’ a? =) E “1 
i- c = Sj) GEN ie im. Co 
| o2,’Vo?, Q ENN 2) exp 07. I, eu (105) 


and substituting (105) into (104), we get 


ja? jae | a” eta aa, 
Pr{R, <R,IS,)= a =e Be Pos ‘I, = ; (106) 


Now, replacing a*, with a*, +a’; +2a,a, cos(9), we get 
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PR, <R,ILS),.8), ,a,,4,0) = 
| a?, fa’. +a°j+2a,a, =) ] f a°. +2a°) +2a,a, 80) 
Q San eee exp} — ———__.— | . 


2? ) ee 
oO’, Gc’, p 20°, 





a, ja°. +a°j +2a,a, cos(8) 
I, o. o* 
(107) 


The dependence on a,, a,, and 8 must be removed by multiplying (107) by the Ricean 


pdfs for the signal and jammer tones, respectively, and by the pdf for 8 and then 


integrating over a,, a,, and 8. As previously mentioned, this must be done numerically 


and can be a computationally intensive operation. In any event, the problem at hand has 


been "solved" in the sense expressed in [7]: 


In many cases we shall consider the solution found if we can express it in 
terms of one or more integrals, even if the integration can only be 
performed by expanding in series (or computing numerically). This is 
called ‘reducing to quadratures', a phrase by which the mathematician 
symbolically washes his hands of the remainder of the task of finding the 
solution. 


This completes the analysis of the receiver performance for independent multitone 
jamming. The final form of the equations that must be evaluated are summarized in the 
appendix, equations B.6 through B.12. Many of these equations are also used for the 


analysis of band multitone jamming, which is addressed in the next chapter. 





VI. FH/MFSK RECEIVER PERFORMANCE WITH BAND MULTITONE 
INTERFERENCE 


Band multitone jamming implies that the multiple, equal power jamming tones 
are distributed randomly across the entire frequency-hopped bandwidth with a specific 
number of jamming tones placed in a frequency-hop bin given that a particular hop bin is 
jammed. For worst case band multitone jamming, the worst probability of bit error is 
obtained when at most a single interfering tone exists at one of the M possible orthogonal 
signaling frequencies in a specific FH band. This section considers the general case of a 
jammer using L jamming tones per hop bin with 1<L<M. We require the probability that 
the jammer places a set of L tones within the signal frequency-hop bin. This probability is 
derived by considering the following numerical example. Suppose that we have a total of 
10 jammer tones(q = 10), and two jamming tones per jammed bin (L = 2). Suppose that 
the number of hop bins is 30(N=30). The jammer tones are grouped into q/L = 5 sets, 
each set composed of two tones. Now, the probability that one of the five sets exists in 


one of the N = 30 hop bins is 5/30. In the general case, the probability that the hop bin is 


jammed is 
< q/L_ sq 
dj=——- = = 108 
Pr(bin jamme ) N “NL (108) 
The probability that the hop bin is not jammed is 
Pr(bin not jammed) = ae (109) 


N-L 


From the total probability law, we find that the overall probability of error is: 


Pr(error)=Pr(L,L > 1)-P.(L,L > 1) + Pr(L = 0)-P,(L =0) 


rai 
= ——_. > -_-——— |. = 
Ki P(L,L21)+|1 isi P.(L = 0) 


(110) 


The probability of symbol error when a hop bin is not jammed, P,(L = 0), is the same as 


the probability of symbol error for noncoherent MFSK and no tone interference 
P.(L=0)<(M-—1)-Pr(R, <R,]). (111) 


This non-jamming case must be considered in order to account for every possible 
jamming scenario. 

The probability of symbol error when a hop bin is jammed is found in (38) to be 
(L21): 


PL) s =-[(L-1)-Pr(R <R,L,S, 8; )+(M-L)-Pr{R, <R 1S, )]+ 
$ ~ M ] Dales? he? jo l 2 I+ iy 
+(1-4).fi-pr{r <R,ILS,,)+(M-L~1)-Pr(R, <R,1)| 
M I a cael 1 2 


(112) 


Hence, combining (110), (111) and (112), we obtain the probability of symbol error 
when (L=1) as 
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le 
— [(L be 1) ' Pr(R, < R,|LS8, 3, \+ (M - IL 5 Pr(R, < R,|1,8; )]+ 





N-L (1-%) [L-Pr(R, <R,ILS, }+(M-L-1)-Pr(R, < Ry) 


G-— | {o - Daou - 5 2) , 





The individual probabilities appearing in (113) are identical to the ones found in the last 
chapter, and are given in (58), (79), (103), and (107). The difference from the 
independent multitone interference is that L is fixed. 

This completes the analysis of the receiver performance for band multitone 
jamming. The final form of the equation that must be evaluated is given in the appendix 
as equation B.13. In the next chapter, numerical results are obtained for receiver 


performance with both band and independent multitone jamming. 
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VII. NUMERICAL RESULTS 


The performance of the FH/MFSK receiver is investigated as a function of the 


; eae F 
ratio of the information signal power to the total interference power, ——. Worst case 


JT 


performance is approximated by selecting q as the largest integer value between 1 and the 


Jr 





ratio . For non-worst-case performance, we select q=100. 


The non-fading, Ricean, or Rayleigh fading channel for the signal and/or the 


jammer tones 1s modeled by appropriately selecting the values of the corresponding 


direct-to-diffuse power ratios, = and m7 . A large ratio value implies little fading, 
a value of 10 implies Ricean fading, and a very small value implies Rayleigh fading. Due 
to programming necessities and in order to avoid divisions by zero, we use a ratio equal 
to 10~ to simulate Rayleigh fading. For purposes of numerical computation, all curves 
are obtained with E,/N, = 13.35 dB and N=1000. The zero thermal noise curve is also 
plotted in each figure to validate the results since this provides an upper bound. For band 


multitone interference with a single interfering tone per hop and zero thermal noise, the 


probability of bit error is equal to [8] 


: ma | 114 
5 ul ers , (114) 


4] 


while for independent multitone interference with zero thermal noise, the probability of 


bit error for FH/BFSK is [9] 


P stale q-! ) tof? -p,) +51) deos" QP. f* 7.) 
ae) IN=V 2\q * 2QNWONeVa 4P, q 4)° 














A. BAND MULTITONE INTERFERENCE 


The performance of FH/BFSK with band multitone interference when the 
information signal is essentially unaffected by channel fading (a ~/20°, = 1000) with 


M=2 and M=8, respectively is shown in Figures 2 and 3. L successively takes the values 


1 and 2 in Figure 2 and the values 1, 4, and 8 in Figure 3. Performance is computed for 





Pp 
both q = 100 and worst case q. When P; 2 P, , or equivalently = < —20dB and q=100, 


* 


Jt 


the worst probability of bit error is obtained for at most a single interfering tone at one of 


the M possible orthogonal signaling frequencies (L=1). For the region where P, 2 P,, the 


worst case performance is obtained when all signal tones are jammed (L=M). Hence, a 
jammer whose power P; is guaranteed to exceed the signal power will do more harm with 
just one tone per bin. However, a jammer with P,; guaranteed to be less then the power of 
a signal tone, will do more harm by jamming all M tones of a bin. For worst case q, the 
worst performance is obtained when L=1 for all power ratios. 

The non-fading signal case for M=2 with all three possibilities of jammer tones 


fading is shown in Figure 4. The results are the same as those in [1]. When P; 2 P, or 





— <-—20dB and q=100, the difference in performance is trivial regardless 


JT 


equivalently 


of whether or not the jammer tones are affected by the fading nature of the channel, even 


42 


when the jammer tones experience Rayleigh fading. On the other hand, when P, 2 P,, or 





Ie 
equivalently = 2 —20dB, the interference tones that experience Rayleigh fading yield 


aT 


the poorest performance. For the worst case q, the difference is trivial for the entire range 
of signal-to-jammer power ratios. Figure 4 is identical to the figure we get if we use the 
exact result equations. In fact, the individual probabilities of bit error used in the exact 
result, as they are described in [1], are identical to those obtained with the union bound 
for the M=2 case. So M=2 will not be treated again by itself since it can convey no more 
information. It 1s important that for M=2 the simpler union bound equations provide the 
exact result. 

The non-fading signal case for M=8 is shown in Figures 5 and 6 for q=100 and 
worst case q respectively. The curves obtained for the exact result and the union bound 
are identical and practically indistinguishable from each other. The previous phenomenon 
with regard to the effect of fading interference tones deduced from Figure 4 is also 
obvious here, but it is even more pronounced, with the non-fading and Rayleigh fading 
curves differing by 5dB. 

The Ricean fading signal case for q=100 and worst case q is shown in Figures 7 
and 8 respectively, where the exact and approximate solution's curves obtained for M=8 


are compared. In Figure 7, we notice that the union bound curves are almost identical to 


Cc 





the exact result ones and only in the region where > —10dB do they appear to have a 


It 
constant difference of 1.38-10~. In Figure 8, the union bound curves behave as expected, 


being practically indistinguishable from the exact result ones for P. < es , while in the 


opposite case they appear to have a slight difference with a maximum value equal to 
1.38-107. The union bound curves again work as an upper bound on the probability of bit 
error. 

Summarizing the results up to this point, since the differences detected are trivial 


ones, for all practical purposes we infer that for M even larger than eight, we can easily 
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resort to the union bound equations to estimate the probability of bit error when the signal 
tones are either not affected at all by the fading channel or the signal tones experience 
Ricean fading. 

The exact result and the union bound curves for a Rayleigh fading channel 
affecting the signal tones and the three possible cases of fading conditions for the jammer 
tones are shown in Figure 9. The difference between the two groups of curves attains a 
maximum value equal to 0.0038. This difference between the exact and union bound 
curves is the greatest noticed up to this point. The union bound curves still act as a bound, 
being higher than the exact result ones, but the accuracy is questionable. This 1s examined 
further in the following figures where we determine if increasing the number of symbols 
M causes a decrease in the probability of error as should normally happen. 

The non-fading and Ricean fading signal curves in conjunction with every 
possible fading situation of the jammer tones is displayed in Figures 10 through 14. The 
value of M is sequentially changed, taking the values 2, 4, and 8. In these figures, no 
anomaly is detected, and the curves behave correctly in all cases with the probability of 
error decreasing as the number of symbols M increases. However, the next three figures, 
Figures 15, 16, and 17, display the Rayleigh fading signal case in conjunction with all 
three possible fading conditions of the jammer tones. In all cases the results obtained are 
incorrect. The union bound still acts as a bound, but the probability of error erroneously 
increases as M increases. 

Summarizing the union bound approach for band multitone interference, we 
see that for the non-fading and the Ricean fading signal case, the union bound approach 
gives very accurate results which are perfect in the first case and deteriorate some in the 
second though still providing acceptable accuracy. However, when the signal tones 
experience Rayleigh fading, the union bound approach should be avoided because the 
accuracy obtained is not sufficient. From a computational point of view, the union bound 
approach is far simpler than the exact approach. The exact approach equations include 


joint pdfs of two random variables expressed with double integrals. On the contrary, the 
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union bound equations are expressed as the sum of individual pdfs of one random 


variable and thus they are expressed in terms of one or more integrals. 


B. INDEPENDENT MULTITONE INTERFERENCE 


In independent multitone interference, the number of tones existing per frequency 
hop bin can vary from zero to M, where it equates to the number of signal tones. As 
mentioned before, the calculation of the probability of error here is time consuming due 
to the inability to further simplify the solution which contains three integrations. 
However, one integration is eliminated for the case where either the signal or the jammer 
tones are not affected by the fading channel, and two integrations can be eliminated for 
the case where both the signal and jammer tones are non-fading. 

The M=2 case for a non-fading signal and various jammer fading cases when we 
have independent multitone interference is shown in Figure 18. As with band multitone 
interference, the results obtained for M=2 are identical with the ones given by the exact 
solution. This is because the exact solution equations, described in [1], again coincide 
with the union bound for M=2, enabling us again to use the union bound with perfect 
accuracy. As in band multitone interference, we notice the performance degradation from 


the non-fading jammer case to the Rayleigh fading jammer case when P, 2 P,. 


The non-fading signal case for M=8 is shown in Figure 19. It is obvious that the 
union bound approach is a valuable tool in this case since the exact result approach 
requires very complicated probabilistic equations with the complexity further increased as 
M increases. However, the union bound approach provides us with a concise expression 
for any value of M. The only means of comparing the curves obtained here with a pre- 
existing result is to get a feeling from the comparable ones obtained for band multitone 
interference case bearing in mind the result obtained in [1], that the probability of error 
for the band multitone interference is poorer than for independent multitone interference. 


In fact, comparing Figures 5 and 19, we see that they have the same shape, differing by a 
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factor of roughly eight. Also, comparing Figures 4 and 18, we see that they differ by a 
factor of roughly equal two. The difference between band and independent multitone 
interference is roughly equal to M. 

Figure 19 is characterized by the phenomenon previously described. For 


P, 2 P; achannel with no fading of the interference tones provides an advantage of up to 


5dB over one with Rayleigh fading. 

The Ricean fading signal case for M=2 with all three possibilities of fading of the 
jammer tones is displayed in Figure 20. As for non-fading signal tones, when q=100 we 
notice a difference in performance for the non-fading and Rayleigh fading jammer tones 


when P, = P, . That difference, however, is now less pronounced. For the worst case q, 


the difference between the respective curves is again trivial. Comparing this figure with 
Figure 3, we again notice that the probability of error 1s less by a factor of roughly two. 

The Ricean fading signal case for M=8 with all three possibilities of fading of the 
jammer tones is displayed in Figure 21. Comparing this figure with Figure 7, we notice 
that the probability of error difference is eight. This confirms the previous observation 
that the decrease in probability of error when we shift from band to independent 
multitone interference is roughly equal to M. 

The Rayleigh fading signal case for M=2 and M=8 with non-fading interference 
for q = 100 and worst case q is shown in Figure 22. The obvious problem is the same as 
the one observed for band multitone interference. The accuracy achieved is not sufficient 
and predicts an increase in the probability of bit error for an increase in the number of 
symbols M. Thus, the union bound approach should be avoided when the signal 
experiences Rayleigh fading. 

The remaining figures, Figures 23, 24, 25, and 26, show the results of 
approximating 9 rather than integrating over all of its possible values. Remember that by 
0 we denote the difference between the phase of the signal and the phase of the jammer 
tones. This is modeled as a uniformly distributed random variable over [0, 27]. As 


mentioned before, in the beginning of this section, for the case of more than one jammer 
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tone per frequency hop bin, we must perform three numerical integrations in order to 
calculate the probability of error: one for the amplitude of the signal tone fading, one for 
the amplitude of the jammer tone fading, and one for the angle 8. To this must be added 
the embedded calculation of the Marcum's Q-function that is also defined in integral 
form. It would be very convenient to find a way to eliminate one numerical integration. In 
fact, the integration over 0 was experimentally found to be the most time consuming one. 


After trying different angle values to get an accurate approximation, it is found that for 


7 
Jo) = 10. m the approximation is so good that the difference between the exact and 


approximated traces can barely be detected. The exact and approximated curves for a 
non-fading signal and various fading of the jammer tones for M=2 and M=8, respectively, 
are Shown in Figures 23 and 24. The exact and approximated curves for a Ricean fading 
signal and various fading of the jammer tones for M=2 and M=8, respectively, are shown 
in Figures 25 and 26. This approximation results in a significant reduction in 


computational complexity and time. 
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VIII. CONCLUSION 


From a probabilistic performance point of view, the results obtained verify those 


of [1] that when P, > P, the effect of fading jammer tones is very small when the signal 
tones are not affected by fading. When P, > P, performance is poorer when the 


interference tones experience fading. When the information signal experiences fading, the 
effect of fading interference tones is less. 

With regard to the union bound approach, it is shown to be very accurate for non- 
fading or Ricean fading signal tones, while for Rayleigh fading signals its use should be 
avoided. 

The application of the union bound for the Ricean fading signal in the case of 
independent multitone interference for M=2 and M=8, that was not treated in [1], shows 
that the probability of error advantage when the jammer uses independent multitone 
interference instead of band multitone interference is equal to a factor of M. Hence, when 
M = 8, independent multitone interference gives probabilities of error eight times less 
than does band multitone interference. 


In case of the independent multitone interference, the calculations involved can be 


7 
much simplified if instead of integrating over all possible values of 8 we let 0 = io 


thus eliminating the most time-consuming integration. 

In summary, the union bound approach taken in this thesis not only simplifies the 
computational complexity involved in obtaining numerical results, but also allows us to 
analytically evaluate situations that are not amenable to solution using an exact approach. 
The principle drawback is the breakdown of the union bound approach when the signal! 


experiences Rayleigh fading. 
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Figure 2: Performance of a FH/BFSK noncoherent receiver with essentially 
no information signal fading for band multitone interference with 
various conditions of fading of the multiple interference tones for 
L=1 and L=2. 
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Figure 3: Performance of a FH/8FSK noncoherent receiver with essentially 
no information signal fading for band multitone interference with 
various conditions of fading of the multiple interference tones for 
L=1, L=4, and L=8. 
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Figure 4: Performance of a FH/BFSK noncoherent receiver with essentially 


no information signal fading for band multitone interference with 
various conditions of fading of the multiple interference tones. The 
results are identical with those obtained using the exact result 


equations. 
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Figure 5: Performance of a FH/8FSK noncoherent receiver with essentially 


no information signal fading for band multitone interference with 
various conditions of fading of the multiple interference tones, g = 
100, using both the exact result and the union bound equations. The 
curves obtained are identical and practically indistinguishable. 
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Figure 6: Performane of a FH/8FSK noncoherent receiver with essentially 
no information signal fading for band multitone interference with 
various conditions of fading of the multiple interference tones, q,,.. 
using both the exact result and the union bound equations. The 
curves obtained are identical and practically indistinguishable. 
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Figure 7: Performance of an 8FSK noncoherent receiver with Ricean fading 


of the information signal for band multitone interference with 
various conditions of fading of the multiple interference tones, 

q = 100, using both the exact result and the union bound equations. 
The union bound curves behave as expected, being a little higher 


than the exact solution ones. For = > -10dB we lose a little bit of 
7ij 

accuracy as the union-bound curves have the same shape with the 

exact result curves but they have a constant difference of 1.3610. 

As in the M=2 case, there is no essential difference between no 

fading, Ricean fading and Rayleigh fading of the jamming tones, 
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Figure 8: Performance of an 8FSK noncoherent receiver with Ricean fading 


of the information signal for band multitone interference with 
various conditions of fading of the multiple interference tones, 
dye. uSing both the exact result and the union bound equations. 
The union bound curves behave as expected, being practically 
indistinguishable from the exact result equations for Pc<Pjy while 


for Pc>PjT they appear a little bit higher, having a maximum 


difference of 1.3410 ”. 
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Figure 9: Performance of an 8FSK noncoherent receiver with Ricean fading 
of the information signal for band multitone interference with 
various conditions of fading of the multiple interference tones, 

q = 100, using both the exact result and the union bound equations. 
Again the fading behaviour of the jammer tone is irrelevant to the 
probability of error. The union bound curves have the same form as 
the ones from the exact result, but appear to display the biggest error 
than in any case treated up to is point, an error having a maximum 
equal to 0.038. 
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Figure 10: Checking the Union Bound accuracy for essentially no fading of the 
information signal for band multitone interference with no fading of 
the jammer tones, for both q = 100 and q,. , for M = 2, 4, and 8. No 
anomaly detected for the P,, using the union bound and P, decreases 
as M increases. 


59 


0.1 


0.01 





1°10 > 
1°10 © 
Pc/PjT(dB) 
o& M=2 
> M=4 
= M=8 


“© M=2 worst case q 
“© M=4 worst case q 
“> M=8 worst case q 


Figure 11:Checking the Union Bound accuracy for essentially no fading of the 
information signal for band multitone interference with Rayleigh 
fading of the jammer tones, for both q = 100 andq,,, , for M = 2, 4, 
and 8. No anomaly detected for the P, using the union bound and 
P,, decreases as M increases. 
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Figure 12: Checking the Union Bound accuracy for Ricean fading of the 
information signal for band multitone interference with essentially 
no fading of the jammer tones, for both q = 100 and q,.. , for M = 2, 
4, and 8. No anomaly detected for the P, using the union bound and 
P,, decreases as M increases. 
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Figure 13: Checking the Union Bound accuracy for Ricean fading of the 
information signal for band multitone interference with Ricean 
fading of the jammer tones, for both q = 100 andq,, , for M = 2, 
4, and 8. No anomaly detected for the P, using the union bound 
and P, decreases as M increases. 
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Figure 14: Checking the Union Bound accuracy for Ricean fading of the 
information signal for band multitone interference with Rayleigh 
fading of the jammer tones, for both gq = 100 and q,,. , for M = 2, 
4, and 8. No anomaly detected for the P, using the union bound 
and P, decreases as M increases. 





30 -20 =i 0 10 
Pc/PjT(dB) 


M=2 
M=8 


M=16 

M=2 worst case q 
M=8 worst case q 
M=16 worst case q 


Figure 15 : Checking the Union Bound accuracy for Rayleigh fading of the 
information signal for band multitone interference with 
essentially no fading of the jammer tones, for both q = 100 and 
que > for M = 2, 8, and 16. Anomaly detected!!! For increasing M 
the probability of error increases. 
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Figure 16 : Checking the Union Bound accuracy for Rayleigh fading of the 
information signal for band multitone interference with 
Ricean fading of the jammer tones, for both q = 100 and q,.. 
for M = 2, 8, and 16. Anomaly detected!!! For increasing M the 
probability of error increases. 
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Figure 17: Checking the Union Bound accuracy for Rayleigh fading of the 
information signal for band multitone interference with 
Rayleigh fading of the jammer tones, for both q = 100 and q,,, , 
for M = 2, 8, and 16. Anomaly detected!!! For increasing M the 
probability of error increases. 


66 


O 
* 


x 
ae 





EDINO=13.39 


29 -20  -15 -10 -9 0 
PcfP|T(dB} 


No fading jammer tone 

Ricean Fading jammer tone 
Rayleigh Fading jammer tone 

No fading jammer tone worst case q 


-- Ricean Fading jammer tone worst case q 
- - Rayleigh Fading jammer tone worst case q 


Performance for zero thermal noise 


4% A 
t. ‘ 
tas, So OOS Bees Petre 
= PEST rele a pet NOY AC en 
== * 


Figure 18: Performance of a FH/BFSK noncoherent receiver with essentially 
no information signal fading for independent multitone interference 
with various conditions of fading of the multiple interference tones. 
The results are identical with the ones obtained using the exact 


result equations. 
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Figure 19: Performance of a FH/8FSK non-coherent receiver with essentially 


no signal fading and various interference fading for independent 
multitone interference. 
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Figure 20: Performance of a FH/BFSK non-coherent receiver with Ricean 
signal fading and various interference fading for independent 
multitone interference. 
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Figure 21: Performance of a FH/8FSK non-coherent receiver with Ricean 
signal fading and various interference fading for independent 


multitone interference. 
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Figure 22: Performance of a FH/BFSK and FH/8FSK receiver for Rayleigh 
fading information signal and non-fading interference, for 
independent multitone interference. Anomaly detected, since 


for an increase in M we notice an increase to the probability of 
error. 
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Figure 23: A comparison of the union bound curves for a FH/BFSK 
noncoherent receiver's performance with essentially no 
information signal fading for independent multitone interference 
with either no fading or Rayleigh fading of the jammer tones, 
with q = 100 and worst case q, obtained by integrating 9 and also 


7 en: 
by approximating 9 with . The results obtained are identical, 


so the approximation can effectively be used simplifying the 
calculations. 
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Figure 24: A comparison of the union bound curves for a FH/8FSK 


noncoherent receiver's performance with essentially no 
information signal fading for independent multitone interference 
with either no fading or Rayleigh fading of the jammer tones, 
with gq = 100 and worst case q, obtained by integrating 9 and also 


7 
by approximating 9 with 10" . The results obtained are identical, 


so the approximation can effectively be used simplifying the 
calculations. 
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Figure 25: A comparison of the union bound curves for a FH/BFSK 
noncoherent receiver's performance with Ricean information 
signal fading for independent multitone interference with either 
no fading or Rayleigh fading of the jammer tones, with q = 100 
and worst case q, obtained by integrating 9 and also by 


i me oe 
approximating 6 with ™ . The results obtained are identical, so 


the approximation can effectively be used simplifying the 
calculations. 
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Figure 26: A comparison of the Union Bound Curves for a FH/8FSK 
noncoherent receiver's performance with Ricean information 
signal fading for independent multitone interference with either 
no fading or Rayleigh fading of the jammer tones, with q = 100 
and worst case q, obtained by integrating 8 and also by 


7 
approximating 9 with 10° . The results obtained are identical, so 


the approximation can effectively be used simplifying the 
calculations. 
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APPENDIX . MATHEMATICAL IDENTITIES AND RELATIONS 


For convenience, all pertinent equations are repeated here. 


A. IDENTITIES 


1. Integrals 


* 1,(2) = exp\ - Jor) J, (jz) =(- 5)" -J, G2) (A.1) 


2 


: I,(z) = —. Jexp|+ Z° cos d — 6) kag (A.2) 
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where S"n(e) is the Laguerre polynomial: 





p 
qn [Bo)_ @ CP a, p> ‘ 
= - =| p! . m-p : = 4 (A.3) 
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|x-exp[-ax ]-1,(Bx)-J,(yx)dx = 5—-exp 7 (2 (A.4) 
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2. The Marcum's Q function and related integrals 








* Q(a,B) = fx- exp - S ) texts (A.5) 
B 
* Q(a,0) =1 (A.6) 
Be Be 
Q(0,0)=exe{ AD 
e Q(a,B) = 1~ (6.2) +ex9{- : )-t6(2.) (A.8) 


° fol zs) R, ( a” oR) I (2 Jap 
= - -ex ——— || = 
Gc C7 S 20°; "Loy 


ai 
O71 +075 a fate. ee] pete 67) +073 a ee o71+072 {pre] 


(A.9) 
B. RELATIONS 
1. Union bound on the probability of bit error 
P, =Pr(R, <R, UR, <R,U...UR, <Ryll) 
< Pr(R, < R,[1I)+ Pr(R, < R,|1)+...4+Pr(R, < Ry ID B.1) 


<(M-1)-Pr(R, <R,]1) 
where by R,, R,,...R,y, we express the random variables representing the outputs of 
branches 1,2,...M, respectively, of an MFSK receiver. 
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2. pdf ofa signal's or jammer's amplitude modeled as a Ricean random 


variable 
2 2. 2 | 2 
f, (a,) = + exp - vie tof 2 -u(a; | (i = c, j) (B.2) 


where u(e) is the unit step function, I,(¢) is the modified Bessel function of the first kind 
and order zero, a”. is the power of the direct signal component, and 20’, is the power of 


the diffuse signal component of the respective tone. 


3. pdf of a signal's or jammer's amplitude modeled as a Rayleigh random 


variable 





2 
a; a ij : : 7 
fy (a, ] ie exp - G2 -u(a; } (i = c, j) (B.3) 
4. Signal-to-noise power ratio of the direct signal components 


2 
Oj : . 
Piao Gi=c,)) (B.4) 


n 





5. Signal-to-noise power ratio of the diffuse signal components 


gat (i=c,j) (B.5) 
Oo 





8] 


6. Pr(R, <R,|1) 





I p 
see 5 B. 
Pr(R, <R, 1) a 2] a (B.6) 
7. Pr(R, <R,{1,Si, ) 
Pr(R Se) = —— ee cow Dei Pipe I A Lite (B.7) 
fe E+E, +2 P E.+€,+2) “lE, +E, +2 . 


8. Pr(R, < R,ILS, ) 


2p, 2p, 
Pr(R, <R, 1,Sj,) | | tee phe 





(B.8) 
eee P. TP; 2 PP; 
a ee |, ee a) Fs AE 
Pa ag LPO ge Gy Bears: 
9, Pr(R, <R,ILS,,5S;,.4..2,,0] 
Pr(R, <R,|1,S}, 5}. ,a,,a;,0] = 
a?, ja’. +a°;+2a,a,cos(6)| 1 a’, +2a°; +2a,a, cos(0) 
Q Go’, 2 o?, ae pias cas 262, © 
(B.9) 


Ae la’, +a’; +2a.a; cos(0) 
0 co, o?, 
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10. Pr(R, <R,ILS,_,S,, ] 


0 Pe a. ra ae +a°j +2a,a; cos(@) | _ 
Gy Gar 


' ¢ + 2a’, pes) 





] 2m 2 “exp 
Pr(R, <R,ILS;,.S;,)=—— | ff 
27 000 (B.10) 
a, ja’. +a°;+2a,a, a°. +a’; + 2a,a; cos(8) 
Ih cos. GC’, 
1 (a,). ,)da, da do 


11. Probability of symbol error for L jamming tones per hop bin, 1<L<M 


P(L)< = [(L- 1)-Pr(R, <R,IL,S, 8, }+(M-L)-Pr(R, <RylLS, )|+ 
(B.11) 
+(1-=) JL-Pr{R, < R,IL,S, )+(M-L-1)-Pr(R, < Rll) 


12. Probability of symbol error for independent multitone interference for L 


jamming tones per hop bin, 0<L<M. 


rp =0-—8 <5) (M—1)Pr(R, <R,|l)+ 


min(M-1,q)[ L- qa ni a 

4 2 Heist [=a ae etl) ae 
Ll gj 

a8 acapRon 


where the probability P,(M) is obtained by letting L=M in P,(L). This yields 
P.(M) =(M -1)-Pr(R, <R,/1,8, ,S; ). 


So 


Ll -sa) POO ne ona | 


c 


min(M-1,q)[ L-1 ae M-1 oe 
st aft 2] 4 


0 a’, Ap eaes + 2a,a, cos(6) 
Cat. ow 


1] a’. +a; +2a.a; cos(6) 
5 EXP a 

















Zz 


a0) > ff ; 


a; ago a i 2a a; cos(6) 4 
I | TE Se Le 
lg o?, 


p, +2, 2/PP; 
soe gah cpol- get) [ 8) 


‘ 2p; 2p, 1+é, P. +P; F 2,/P-P; 
1) CU fave. +8, f2+k, +e} 248, +8, M248, +8) 248, +8, 


|= 





























0 a’; a’, +a; +2a.a; cos(@) 
orn Gin 


= ir dex a £.(a,)-£, (a,)da,da a6 
000 


p - 
; a, fa’. +a’; +2a,a, cos(8) 
Ol a2. 


(B.12) 


13. Probability of symbol error for band multitone interference for L 
jamming tones per hop bin, L21 


L 
a [leone cns,s,)otuev-n(e <x 
- e 
~N-L ( 
1 





Ps(L T 
-) |L-P:(R, <R,ILS,, )+(M-L-1)-Pr(R, <R,I) 


1-54] {o-9 oe) exg(- 2 a.) 


SO 








q e 
N-M 


a dae +a’; +2a,a,cos(@) 
Q > 9 2 _ 
Oo On 


a +a’; sce) 
2 — 7 





Ps(L > 1) < 





a 


l 
i 270 0] — 5 
sili E a*. ta’; +2a.a; 248) ay 
I —a ee 
: So”, or 


\ f, (2,)-f,, (a; }da,da 40 J 
no leank ta 
ee ej +€,+2 °°? Sj re _ es +2 
VAP i Os oe 
Q 2+EL +8, eae —24+é, +E, 
E 
q oye 2 P.P J B.13 
+i =) exo(-=2 Tee, fa AE een be | | ie 
a oe re OP ra) 


+(1-585) {o-n; ca - a 
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14. Relation between probability of bit error and probability of symbol error 


M/2 


P. =——— 
7M = 1 


-P., where M is the number of symbols (B.14) 
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